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Abstract 



The mathematical description of laboratory fusion plasmas produced in Tokamaks is still challenging. Com- 
plete models for electrons and ions, as Vlasov-Maxwell systems, are computationally too expensive because they 
take into account all details and scales of magneto-hydrodynamics. In particular, for most of the relevant stud- 
ies, the mass electron is negligible and the velocity of material waves is much smaller than the speed of light. 
Therefore it is useful to understand simplified models. Here we propose and study one of those which keeps 
both the complexity of the Vlasov equation for ions and the Hall effect in Maxwell's equation. Based on energy 
dissipation, a fundamental physical property, we show that the model is nonlinear stable and consequently prove 
existence. 

Key words: Vlasov equations; Plasma physics; Kinetic averaging lemma; Maxwell's equations. 
Mathematics Subject Classification: 35B35, 35L60, 82D10 

1 Introduction 

To describe the behavior of ions population in hot plasmas, it is very classical to address, at least at theoretical 
level, a Vlasov equation coupled to a non-linear Poisson equation which defines the electrostatic field. When the 
magnetic field B is an external datum, for the ion distribution function /(<,a;, v) (at position x and velocity v), 
one addresses the following system 



where rie = ne{t, x) is the electron density. To is the mean electron temperature assumed to be constant and A is the 
Debye length (a characteristic constant of the plasma). This system is a classical one, indeed the electrostatic field 
may be approximated by — To V In rie (at least when the electron temperature is constant) and the Poisson equation 
is nothing but the Gauss relation applied to this field (cf. [7] for example) 

The mathematical understanding of this kind of kinetic system has made important progresses with the proof of 
existence of global weak solutions in the large of the Vlasov-Maxwell system in [TT] and of the Vlasov-Poisson system 
P^I25j . However, the mathematical description of laboratory fusion plasmas (such as those produced in Tokamaks) 
with this kind of ion kinetic models is still a major challenge, in particular to deal with a time scale compatible with 
the evolution equation for the magnetic field. As a matter of fact, complete systems of Vlasov-Maxwell type which 
account for all scales of electro-dynamics are not relevant since the velocity of ion waves is much smaller than the 
speed of light: it is well known that at the time scale of the ion population it is convenient to neglect the current 
of displacement in the Maxwell equations (as in magneto- hydrodynamics models cf. |201 [21]), that is to say the 
electric current J is assumed to satisfy 
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J = V AB. 



The aim of this work is precisely to propose, justify and study a model which couples a kinetic equation for the 
ions and an evolution equation for the magnetic field (as those used in magneto- hydrodynamics). The unknowns 
are the ion particle density /(t,.T,v), the magnetic field B(i, a;) and the electron density ne{t,x) and they satisfy 



— A'^A In rie = n/ — Tie, 

dB ^ / 1 

— - V A — Ti/uj A B 

ot \ne 

ot ov 



VA I — J AB 

rip 



To J — n,Ur _ _ , 

- — Vne+ — A B + V A B / 



V A (77V AB) = 0, 
0, 



[ VB = 0, 
where the following notations are used 

ni{t, x) - 



(a) 

(b) 

(c) 
id) 



(1) 



f{t,x,v)dv 



(the number density in ions). 



(2) 



nj{t, x)uj{t, x) = / /(t, .T, v)v(iv (the macroscopic velocity of ions). 

and r] denotes a strictly positive bounded function corresponding to the plasma resistivity. 
Notice firstly that the electric field E is given by the relation 

rieE = -TcVne - n/U/ A B + J A B + 77V A B 

which is one of the classical forms of the generalized Ohm law [71 [31 [H |T3] (it includes the term J A B related to 
the Hall effect); so, (II|)-b) is exactly the Faraday equation 



(3) 



(4) 



dB 

at 



VAE==0. 



We have chosen here a rescaling, such that the ion mass is set to 1 and the electron and ion charge are also set 
to 1. Moreover the scaling of the ion and electron density is such that the characteristic value of n/ is equal to 1. So 
with this scaling, the Debye length A is defined by A'^ = TqEq where £0 is vacuum dielectric constant and equation 
([T|)-a) reads also as 

— eO^cAlnTle = 71/ — Tie 

which is Gauss relation applied to — TeVlnrje, the dominant term in the electric field. In the rest of the paper, we 
use the following notations 

E°(t,a;) = -Tc— + — (J -ri.uj) AB and F(t, cc, v) = E° + v A B. (5) 

Tie 77e 

Here F is the kinetic force field which appears in equation (HJ-c). 

Remark 1. It is fundamental to notice that the Debye length is a small quantity in many situations; it may be 
proved that in the limit A — > the solution Tie to the non-linear Poisson equation satisfies ?ie — s- ti/ and many studies 
have analyzed this singular limit in different contexts 151 [T5|. Then, if we make the approximation Tie — nj, the 
Ohm law (|4]) leads to the classical formula 



E + uj A B 



-TcV In Ti/ 



1 

Tl/ 



-(V AB) AB + T/V AB 



and from ([T]), we recover formerly the following system 

nr> / 1 

— - V A (u/ A B) + V A — (V A B) A B) 

ot \ Tl/ 

dt dv 



- V A (lyV A B) 
[(E""^ + V a B)/] 



- 0, 



0, 



n/(E" 



u/AB) = -TcVti/ + (V A B) A B. 



Despite its apparently simpler form than the system ([T]), it is not used in physical literature, up to our knowledge. 
As a matter of fact even if the magnetic field is neglected, the equation 



df 



d 



^ + ^-^Z + 7^ [(E'"" + ^ ^ B)/] = 0, Tl/E'"" = -TeVTl/ 
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is probably ill-posed from a mathematical point of view (nevertheless |18j has noticed that using some arguments of 
the paper |22| , one may find a weak solution on a small time interval) . Moreover it would be hard to give a meaning 
of the quantity u/ = J^^ /(v)v(iv/ J^^ /(v)(iv which appears in the magnetic equation (it is hardly possible to 
derive a bound of nj away from below). Notice that this kind of model has been addressed in jTB] and jl] but 
the framework is a little different since a velocity diffusion term (of Fokker-Planck type) is accountered for in the 
kinetic equation. 

Back to |T]), from a mathematical point of view, one can show that Ue is naturally bounded away from zero due 
to the nonlinear elliptic equation ([ij-a and the quantity ^^uj = Jjjj3/(v)v(iv may be defined in some Lebesgue 
space. 

Our study is aiming at the analysis of model ((l} in a bounded domain; therefore we have to specify the boundary 
and initial conditions. 

Boundary conditions Let fl he a. bounded domain of M!^ of class C^'^. We consider the model ([T]) for v G M"^, 
X £ Q and t > 0, with the following boundary conditions 

nj. • Vnf.{t, x) ~ X £ dfl, (a) 

AB(t,x) =0 xedn, (b) (6) 

/(x,v-2(vn,)n,) = ./(2:,v), v £ R\ x £ dfl. (c) 

The zero flux boundary conditions ©-a) on electrons enforces in ([ij-a) the global neutrality of the plasma 

Tiedx = / riidx. 
n Jn 

The condition (b) on the magnetic field is the simplest one and a more realistic condition is discussed in section 
Our result will be still true with more realistic boundary condition, of impedance type for example, or also a non 
homogeneous boundary condition such as 

n^ AB(t, x) = n^ ABimp, x £ dVl, 

which is much more relevant in the context of confined plasmas in Tokamaks. Condition (l6])-c) is the so-called 
specular reflection; it implies the no-slip boundary condition 

u/(t,a;) -n:^ = 0, x £ dn. (7) 

Our main result is the following. 

Theorem 1 (Nonlinear stability). Assume (|14p and that the initial datum satisfy (|f f p - ()f 3p . From any hounded 
family of solutions /'^,B'^,n| to ([T]), ([6|) we can extract a subsequence that converges to a weak solution with finite 
energy. Moreover, we have for all T > 0, 

B^ ^ B e L°°(0,r;L2(r2)) nL2(0,T;L«(r2)) strongly, q < 6, 

t f £ L°°(0,r;Li nL°°(r2 X R^))) weakly, 
J^ = VAB^^J \nl?{{Q,T)y.VL) weakly, 
nl — >■ He in L'((0,T) x il), strongly and for some constants < if_ < n^. < -R'^-, 

(L'?(0,r X f7))3 strongly, 

for all q such that f < g < 5/4. 

The first tool towards a mathematical analysis is to state the energy balance for the full system; this is performed 
in section[21 It indicates that the magnetic resistivity in equation (HJ-a) is useful (if not necessary) so as to control the 
current J and subsequently the term J AB describing to the Hall effect. Also the classical tools for elliptic equation 
with Neumann boundary conditions will lead to control rig away from below. But more technical ingredients (namely 
a kinetic averaging lemma in a L"^ space) are used in the proof in order to pass to the limit in the nonlinearities; 
this is the purpose of section [3] and of several appendices. We come back on more realistic boundary conditions 
in section U] so as to include confinement. We also detail a constructive splitting procedure which can be used to 
design an approximate solution with the same energy law. 
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2 Preliminaries 

The first and major ingredient in our approach to nonlinear stabihty is the energy balance. We state it here together 
with considerations on the physics sustaining the model. 

2.1 The energy balance 

We now turn to the study of energy dissipation for the system ([T]) with homogeneous boundary conditions ([6]) . We 
introduce the ions kinetic energy and the magnetic energy 

as well as the free energy / {ne{t) hi rig (t) — ne{t) + l)dx and the total energy which reads as 

Jn 



Stotit) = Slit) + Srnit) + ^ [ \V\nneit)\^ dx+ [ {neit)\nneit)-ne{t) + l)dx. 

^ Jn Jn 

f 2 

Notice that the integral — / |Vlnne(t)| dx corresponds to the electrostatic energy. We now establish the 

J o 

Proposition 1 (Energy dissipation). Classical solutions to (Jlj, ([6|) satisfy the energy dissipation relation 



(8) 



d_ 

dt 



Si + £m + -z- I |V (Inrie)!^ da; + / [nelane — + 1) dx 





V AB 


Jn 





dx. 



(9) 



Proof. Notice that, as usual, the free divergence condition ([I])-(d) is not needed for the energy identity. We first 
consider the ion kinetic energy. Using the specular reflection condition, we compute 



dt 



n Jr^ 



F(t, X, v) • vdv 



n JR3 



— (VAB)AB + (v-— u/)AB-TeV(lnn,) 1/ 

Hp n„ 



■ -vdvdx 



= ^ A B - n/TeV (Inne)^ • uidx. 

Now, we turn to the magnetic energy and first recall the definition of E in Thanks to a classical tensorial 
identity, we find 

d f c^B f f 

—£m = / ^ ■ Bdx = - / V A E • Bdx = - / [E-VAB + V-(EA B)]dx 
at Jn ot Jn Jn 



ni 



u/ A B + 7/V A B 



(V A B) dx ~ n.^-(EA B)dx. 
Jon 



According to the boundary condition ©-(b), we have n^j • (E A B) ~ — E • (rixAB) = 0, then the boundary integral 
is zero. Now, thanks to identity J • (u/ A B) = uj • (J A B), we get 



^fm = - / — u/ • (J A B)dx - I rj 

at Jq_ Tie 



V AB 



dx. 



At this stage, we have obtained 





V AB 


In 





dx 



J^niV (Innp) -Uidx - J^T] 



V AB 



(10) 



dx. 
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To continue, we recall the conservation of mass on rt/ 

dn 



dt 



- + V ■ (n/u/) = 0. 



Then, we may use the boundary condition ([7]) to obtain 

— / n/V (Inrie) • U/da; = / (Inne)V • (n/U/) 
Jn Jn 

= - / (Inne)— -^dx 
(In rig) 



- A^ — A (Inne 



dx 



d 
'dt 



(rig (In Tie) ~ ne)dx 









/ 


dt 


JO 



|V(ln7T,e 



-dx 



where we have used equation ©-(a) once differentiated in time. 

Altogether, we obtain the relation ([9]). □ 
2.2 Physical considerations 

Let first focus on the momentum balance. To find it, we multiply the kinetic equation ((ij-c) by v, integrate over 
R'^ and define as usual [6] the ion pressure tensor P/ through 



y V ■ (v (g) v/(v))dv = V • (n/U/ (g) u/) + V ■ P/. 



Then, we get 

dt 

that is to say 

a(?i/U/) 



dt 



V • (n/u/ (8) u/) + V • P/ + Te -m + — n,u/ A B A B - n^ui A B = 0, 

Tie 



V • (?i/u/ u/) + V • P/ + ToVTie - J A B -^-^^ — [TeVne + (?i/U/ - J) A B] 



On the left hand side one recognizes the classical momentum equation which appears in the MHD modeling with 
the magnetic pressure tensor. Let us stress that with the electric field E° introduced in ([5]) and according to ((Ij-a), 
the r.h.s. term reads as follows 



V • E" + V • 



niuj — J 



AB 



E"=eoV- -|E°|2-E"®E° 



A^E^V 



71/ Uj — J 



AB 



Thus we may see that this equation is in a conservative form only when the Debye length vanishes. This is due 
to the asymptotic expansion motivating ( [1]) ; we have chosen to keep energy dissipation and smoothness for rig to 
the expense of momentum balance. In order to ensure exact conservation momentum, a possible route is to use a 
modified Poisson equation as follows 



-A-'Alni 



nj ~ He — A^V ■ 



AB 



which corresponds to an approximation at the same order but which looses nice properties on Tie- 

Remark 2 (Electron momentum balance). We see that E" is equal to the electric field E up to the resistive part 
77V A B. Let us stress that (V A B)— ti/U/ = J — njUj is the electron contribution to the electric current and @ 
may be interpreted as the electron momentum balance equation when using the massless electron approximation. 
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Remark 3 (Friction). In order to account for the friction between ions and electrons in the kinetic equation, we 
could add a friction term to equation (HJ-c) and arrive to 

df d d 

(recall that the coefficient rj is related to the collision frequency of electrons against ions). Also, for the sake of 
compatibility, instead of equation (HJ-a) , the Poisson equation would read 

- Alnne + 7^V(77J)) = n/ - n^. 

^ C 

3 Proof of the stability theorem 

We give a family of initial conditions that satisfy 

/^(O, x, v) = r"'^(a;, v) > 0, /™'^ is bounded in n L°°(l] x R^), (11) 

B-'(0,x) =B"-(x), V^-W^' = Q. (12) 

From the Lemma [1] in Appendix VK\ and the explanations in Appendix IB] we see that if 5/(0) is uniformly bounded 

then the free energy at initial time / jV^: (In 71^(0, x))|^ (i.T + / n1[0^x)\rLn1{Q,x)dx is also uniformly bounded. 

Jn Jn 

Then we only need to assume that 

sup [£iiO)+£„-,iO)] < oo (13) 

0<£<1 

to have a total energy uniformly bounded for all time. Moreover we assume a control on resistivity as 

0<77mi„ <r;eL°°(r!). (14) 

Associated with such initial datum, we consider a sequence of strong solutions , B*^, and of ([T|) with boundary 
conditions ([5]). 

3.1 A priori bounds 

Before we begin the proof of Theorem [1] we recall several general a priori estimates that are used. The energy 
bound (|13p is of course at the heart of our analysis because it allows us to deduce directly from the estimate (j9]) the 

Proposition 2 (Bounds derived from energy). Under the assumptions (|lip - (jl4p . the sequences B*^, and 
satisfy 



sup sup / / /f (t,a;,v)|vp(ivda; < c», (15) 

o<e<i te[o,oo] Jr3 Jn 

sup sup / \B''(t,x)\'^dx < oo, (16) 
o<e<i te[o,oo] Jn 

sup sup / (\W^{\n\nl\)\^ + nl\n\nl\+nl) dx < oo. (17) 

0<e<l t6[0,oo] Jn ^ ' 

Indeed, for (|17p . we use the fact that the map a; G M i— > x\i\{x) — x is bound from below and thus all the terms 
in 19]) are bounded by the initial energy plus a constant. 

This provides a control on the kinetic energy for the ion distribution that is essential to obtain further a priori 
bounds 

Proposition 3 (Estimates for the kinetic densities). For strong solutions one has e L°° (O, oo; nL°°(il x M?)') 
and 

\\r{t.:-)\\p = \\n'{-,-)\\v vi>o, vpe[i,oo], 

||nf(t,-)ll5/3<qi.r'^(-,-)ll^' (^I^I^Jit,x,^)H'dvdxy\ 

||nfuf(i,.)||5/4<C||r"'n-,-)ll^'f/ / fit,x,v)\v\'dvdx)^\ (18) 

\Jn JR3 / 
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The bounds of follow immediately from the observation that divv-F = 0. The bounds on n'j and n'juj follow 
from interpolation inequalities that we recall in the Appendix [X] 

Next we use the L^^^ integrability of ni to obtain 

Proposition 4 (Estimates on the electron density). The electron density satisfies for some constants > 0, 
> (depending on the bounds stated as now) 

< K_ <ne{t,x) < K+. (19) 

This is an easy consequence of elliptic regularity and we give a proof in the Appendix |Bl it is just a combination 
of Lemma [1] and Lemma [51 



We now come to the magnetic field. We have the following Proposition : 
Proposition 5 (Estimate on the magnetic field) . The family of magnetic fields satisfies for a uniform constant C 



sup 







JO Jo. 





2 

dxdt < C, 



Proof. The control in L^(0, oo; Hcuii(f^)), follows from the energy dissipation in Proposition [TJ Then thanks to 
Theorem[21 the family is bounded in L'^(0, oo; H^(ri)'^), and then, thanks to the Sobolev injections, the family 
B"^ is bounded in L^(0, oo; L^(ri)'^) . We refer to Appendix [Cl for precise statements and references. □ 

Our last task is to show that we have enough bounds to define the electric and force fields. The four terms that 
compose the electric field 

W = -TeV In 72l - A B^ + + 77V A B^ 

are respectively uniformly bounded in the spaces 



TeVlnTi^ e L'^{L^) (energy inequality), 



AB" e (L^°/^^); this follows from (HH), (HH) and ^ and the fact that, according to Holder inequality. 



we get the bound ^ [\\AB\\.,,/,,]'dt < ^ [\\A\\^/, \\B\\,?dt < sup \\A\\l/, ^ [\\B\\,]'dt. 

je A B"^ 

• e L( (ii^/^) nL?(ii,) from the Holder inequality with (|19p and the a priori estimates in Proposition [5l 

n| 

By interpolation, we also find that 

nl p S 6r 

In particular ^^^^ & Lf /"(if/^^). 

• 7]\7 A B" G 2; (energy dissipation). 

Therefore it is well defined and the four terms will have a weak limit in Lebesgue spaces. 



For the kinetic force field, 

AB^ 



F''(t,x,v) -ToVln?!^ + v - -^-m AB 

V "I / 

the same integrability. locally uniformly in v, holds true because these are the same terms and thus they are bounded 
as for the electric field at the exception of v A B^ which is well defined thanks to ([T5)) . 
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3.2 Space-time compactness 

The bounds mentioned before allow us to extract subsequences that converge weakly and our task is now to prove 
several strong convergence results. 

We first observe that, from Proposition [5] and (fT6|). we may extract a subsequence, still denoted by such that 
— ^B strongly in LP(0,T;L«(17)^), 1 <p < 2, 1 < q < 6, and 1 <p < oo, 1 < g < 2. (21) 

e— s-O 

This is because (see Appendix [C|) thanks to our a priori estimates and Theorem [21 the field {B^} is bounded in 
Li(0,T,Hi(f])3), and then, thanks to Rellich Theorem, {B^(t, •)} is relatively compact in h^{nf for all t e]0,T[). 
Compactness in time then follows from the Simon-Lions- Aubin lemma (see [28]) because, as proved at the end of 
section p.ip . we control L[(LP) integrability of the electric field and 

dt 

According to the kinetic averaging lemma recalled in Appendix[Dl we know that after extraction of a subsequence, 
for every test function = ipi'^) compactly supported, there exists a weak limit /* of such that 

ii^n^ii^n a.e. and (Vv/^) ^ (^.v/*) a.e. 

Now, since / / | (i-'f^) f^^dxdt and / / | (ipv-f^) f^'^dxdt are uniformly bounded, according to classical inter- 
Jo Jn Jo Jn 
polation arguments we have for given exponents q and r (with 1 < g < 5/3 and 1 < ?' < 5/4) 

/ I ii^n - (V'/*) l^dxdt -> 0, and [ [ \ (^'v/") - (^Av/*) ^dxdt 0. 
Jn Jo Jn 

Therefore we get 

n/, strongly in L«(0, T; Q), (22) 
n/uf — > njui strongly in L'"(0, T; n)^. (23) 
Next for the electronic density we may extract a subsequence, still denoted nl such that, 



Tie a.e. and strongly in L«(0,r; iiy. (24) 



This is a consequence of the bounds (|17p and ([T9|) and the strong convergence ([22)) . 

3.3 Passing to the limit 

We recall that a weak (distributional) solution to ([!]) is defined by testing against smooth test functions ^'(t, x), 
^{t,x,v) respectively for ([ij-b) and ([Ij-c) and satisfying the corresponding boundary conditions as in ([6]). 
After integration by parts, we obtain respectively for lU-b) and Q-c), the definitions 



(0,T)xf2 



dt 



+ E''{t,x) -"^ A'iit,x)\dtdx = [ B"^^(x, v)*(0, x), (25) 

Jn 

, V • V^^{t, .T, v) + F=(i, X, v) • Vv<E>(t, X, v) f^{t, X, ■v)dtdxdv 

(26) 

= //' r"'-(x,v)<i>(0,a;,v). 

The elliptic equation ([IJ-a) is more standard and does not yield difficulties, hence we do not consider it here. In 
order to conclude the proof of Theorem [1] our purpose is to show that the limit as e — > of the various unknowns 
B^, W{t,x), F' and still satisfy these equalities. 



and 

-9$(i,a;,v) 

(0,T)xOxK3 
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Additionally to the strong convergence results stated in (|2ip -(|M |) . we are know other weak convergences for 
functions of interest. Firstly, we have 

This is enough to pass to the limit weakly in the electric fields because the nonlinear terms are always formed 
of either strongly convergent terms or multiplied by a term that converges strongly. As a conclusion, we have 

— ^ E in L2"/" (0, T; L30/29(^^))^ 
The same applies to the force field. However for later purposes, it is better to use the splitting 

and to notice that, as before, 

go.e ^ .^^ L^"/" (0, T; L3°/29(r!))' 

while the other term converges strongly. 

These observations allow us to pass to the weak limit in equations ([Ij-a) and ([Ij-b), on the electronic density 
and on the magnetic field. 

For passing to the weak limit in the third equation (UJ-c), let us introduce two smooth test functions x ^-nd if 
with compact support respectively in [0,T) x fl and in R'^ ; then we have to deal with two terms. The first one 
which reads as ^ 

/ / / B''-' Av ■ ^{v)x{t,x)r(t,x,v)dtdxdv, 
Jo JnJm 

can be treated as usual because the magnetic field converges strongly as stated earlier. The second one reads as 
follows ^ ^ 

Jo JqJm=> ov Jo 

with ipi = ; we need to use the kinetic averaging lemma which is stated in appendix [Dl It states that, after 
extraction, velocity averages converge strongly 

xm'){t,x)^x{^^n ill L«([0,T]xr!) 

for a given value of the index q, with 1 < q < p. Therefore they converge almost everywhere and this is enough to 
pass to the weak-strong limit in the term J^^ xE*''^ • (i^f^) dtdx. 

This concludes the proof of the Theorem [TJ 



4 Non homogeneous magnetic boundary condition 

In view of applications to the physics of confined plasmas, it is worthwhile considering, instead of ©-(b), a non 
homogeneous boundary condition for the magnetic field 

n:,AB(<,a-) ABin,p, x e dfl, (27) 

The imposed magnetic field Bi^p{x) is given and can be a function of time as well eventhough we restrict our 
analysis to space dependency only. In order to introduce this boundary condition in the problem, we assume that 
Binip can be defined globally as a smooth function in the closure of O. It turns out that a convenient regularity 
assumption is 

||Bii„p||i,oo < C. (28) 



9 



We will write B = Bi„ip + Bpcrt where Bport is the perturbation. Using these notations the system (IT])-(a,b,c) is 
rewritten as 



— A^Aln7ie = nj — rig, 



(a) 



dB 



pert 



dt 



7^ + v.V/ + — 
ot ov 



V A ( —m Uj A (Bimp + Bpcrt) ) + V A ( — (Jimp + Jport) A (Bi^p + Bport) 



-V A (7/(Jinip+ Jpcrt)) = 0, (6) 



(_Z£Vn,+ ^^'-P + -^P'-'-^^ "-"^ A (B,„p + Bpert)) + v A (B,„p + Bp^,) ) / 



= 0, (c). 
(29) 



where the total current is J = Jimp + Jpcrt with 



Jimp — V A Bimp and Jpcrt — V A Bport- 



We now use the boundary conditions 



■ Vne{t, a;) = 
AB(t,a;)pcrt = 



X G dil, (a) 
X G dn, (b) 

(c) 



(30) 



/(x,v-2(vn,)n,) =/(x,v), v G M^ x e dil. 
Let us define the perturbed magnetic energy and the total perturbed energy 

£r'it)^l I iBport(i,x)pdx, 

f Pf ^Ei + £P"* + ^ / |V. (Inne)l' dx+ [ {rip Inn, ~np + l)dx> 0. 
^ Jn Jn 

The energy balance is modified by the imposed magnetic field and we have 
Proposition 6. Classical solutions to Il29\) - I[3d\} satisfy the perturbed energy dissipation relation 



d ppert 



/ V 


Jpcrt 


Jn 





dx + S 



where the source is 



^1 Jimp ' Jpcrt^^ 



"Jimp A (Bjmp ~t~ Bpcrt) ' Jpcrt^^ 



1 



Jimp A (Bimp + Bport) • "7 VLldx. 



Proof. Performing the same manipulations as in the proof of the energy identity ([9|) , we obtain 

= (Jimp + Jpcrt) A (Bimp + Bpcrt) - n/ToV (luTle)^ • Ujdx. 

Taking the scalar product of the magnetic equation against Bport and using the homogeneous boundary condition 
(1201)- (b), we obtain 

-rSm"*^ — ~ / — U/ • (Jpcrt A (Bimp + Bpgrt)) dx — [ (Jimp + Jpert) A (Bimp + Bpgrt) ' Jpert^a; 

dt Jn Up Jn Up 

/ V (Jimp ^" Jpcrt) ' Jpcrt^^ 

= — I — -U/ ■ (Jpcrt A (Bimp + Bpcrt)) dx ~ [ Jimp + A (Bimp + Bport) " Jpcrt^a; — / ?7 (Jimp + Jpcrt) ' JpcrtC^a^ 

Jn Jn Jr3 



Therefore one gets 



dt 



pcrt\ 
tot / 



TcV (In rie) • {niUi)dx 



' pert 



dx + S 



which is very similar to PU)) . The rest of the proof is unchanged. 



□ 
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Proposition 7. There exists a constant K > depending only on the initial data and on the constant in I128\) such 
that the perturbed energy is bounded for all time t < T* ^ log 11' ^ 



llJ.mpllc 



Remark 4. Our analysis of weak stability can easily be extended to this non- homogeneous boundary condition 
for t < T*. This estimate expresses the interest of a good control on the imposed current. Indeed the smaller is 
||Jimp||oo, the greater T* is. 

Proof, (of proposition [7]) • We observe that the new terms in S are proportional to Jimp. Therefore two cases occur. 

First case: Jimp = 0. This idealized case might be encountered in Tokamaks: for example if Bimp = FV0 
where 9 is the toroidal angle and is a constant. Physically it corresponds to an imposed exterior magnetic with 
vanishing current [TU]. The magnetic lines of Bimp form a ring. In this case S" = 0, so the perturbed energy 
dissipation relation has the same form as It turns out that the perturbed energy is bounded for all times, which 
indeed corresponds to the claim since Jimp — 0. 

Second case: Jimp 7^ 0. This situation is more realistic in Tokamaks, since the magnetic lines form an helix 
or at least must be close to helicoidal geometry [TO] . This is the general case. 
Let tr > be a given positive number. The energy identity writes also 



/ V 


J pert 


In 





^dx + S. (31) 



Our goal is to control the source term S* in i? as much as possible by —aE^^l'^ — fy||Jpcrt||2 and to control the 
remaining part with a Gronwall technique. 

Using the assumption (pS)) . the fact that r] is constant and the Holder inequality between the conjugated spaces 
L^(ri) and ^^(Jl), we can write 



15"! < (ai||Jport||2 + Q'2||Jpcrti|2||Bpcrti|2 + ^2 1 1 Bimp || 5 || U/||5 + Q!2||B 

port II 5 11^/ U-/||5 J 1 1 Jimp 1 1 oo 

with ai ~ T] and Q'2 = ||?T'^^||oo- Let e > be an arbitrary real number. Then 

M T II II ■¥ II lln II '-^1 '-^l ^2 II T ||2 ^-^2 ||T-> Il2 

"iPportlb + a2||Jpcrt||2||Bport||2 < ^ + « ^ l|Jport||2 + 1 1 Bpcrt || 2 ■ 

The identity ([M]). written as U/ljs < a^Sf , implies successively the controls 

4 

IIBimpllsll"/ u/ll 5 < asSf 

llBpertllsl!"/ u,||| < a, (||Bp„.t||2 + ||Jpc..t||2)^/ < ae {S^'Y £] + ||Jpcrt||^ + 
So the right hand side in (|5l]) is bounded by 

for some constants /3i,2,3- Since the electronic density is bounded and < + fP™* < ^t^t'^ by construction, we 
also have 

for some constants 71,2,3 which do not depend on time. Since < 1 + for positive y, we get the more compact 
form 



13 

tot 



for some constants Si,2 which do not depend on time. 

We see that the right hand side is more than linear with respect to £f^^^ due to the power . As a consequence 
this inequality cannot prove that e~'^*£^^^^ or £^^^*' is bounded for all time. 

Next we wish to obtain a evaluation of the time of existence with respect to || Jimp||oo- We set u{t) = e-'"*£:P"f . 
One can simplify the inequality as 

U'{t) < (^Si+S2UT^^ IIJimpllooeT^*. 
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Rescaling of the time variable as dr = ||Jimp|iooe that is 

10, 



■||Jimp||oo (eiS* 



yields the inequality -^u < Si + S2Ui<> . It is finally convenient to define v = Si + S2ui<> so that 

13 3 3 16 

v'{t) = S2—uiou'{t) which yields v'{t) < S^vt^v = S^vt^ . 

Therefore —■§pv^'^ S4 ~ j^^s- which implies v^^{0) — v^'^{t) < S4T. It yields 

?;i3(t) < —3- 

V 13 (0) — S4T 

3 

which is valid for r < t* = Going back to the time variable t, the solution is defined for t < T* where 

10 



yPimplloo [e^o- 



The proof is complete. 



□ 



5 Construction of an approximate solution 

In order to complete our theory, we now detail how to use the so-called splitting strategy, which is a constructive 
method, for the design of an approximate solution to the system 

' — A^Alnjie = ni — n^, 



- VA ( :^ni u/ A B ) + V A (;^J A B)+V A (77V A B) = 0, 

= 0, 



T„ J — n.Ur „. „ 1 „ 

( ^ Vne+ — A B) + V A B / 



(a) 
(b) 
(c). 



(32) 



The idea is clearly inspired from numerical methods. It consists of a convenient splitting strategy a la Strang, 
together with the linearization and freezing of certain coefficients a la Temam. The main point is to decompose 
the total system in simpler parts which are conceptually easier to solve or easier to analyse, preserving at the same 
time the decay of the energy identity 

irie)! dx + (rie In Tie — TT-e + 1) dx. 

In 

Let At > be a time step which is ultimately destinated to tend to zero. We consider that 

/(ifc) and B{tk) 

are known at the beginning of the time step t^. = kAt. We restrict the presentation to the core of the method. This 
constructive method also provides additional insights into the mathematical structure of the model. 



£tot=^l f f f{t,x,v)\v\^dwdx+l [ \B{t,x)\^dx + ^ [ |V,(ln7 
^ Jn ^ Jn ^ Jn 



5.1 Vlasov-Poisson 

One first solves during the time step At 



-A^Aln Ue 

^ = 
dt ' 



9/ 

dt 



v.V/ 



d_ 
dv 



T 

rip. 



= 0. 



(a) 

(b) 

(c) 



This is a non linear Vlasov-Poisson equation which can be considered as standard even if we know very little 
mathematical literature about it. It is easy to show that regular solutions preserve the energy. This procedure 
defines a new solution 

f*itk + At) and B*{tk + At) =B{tk)- 
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5.2 Magnetic part, first stage 

For convenience we split the magnetic part of the equations in two stages, the first one which is fundamental, and 
the second which deals with less involved terms. The first stage writes 



dB 

'dt 

dt dv 



VA ^711 ui A Bf^.o,gjJ + V A (;^J A Bf^„^g^)+V A {fj\7 A B) = 0, 



ib) 



— ^ Bfj.02en) / 



= 0, 



The initial data is provided by the previous step of the algorithm 

f°{tk) = + At) et B°(tfc) = B*{tk + At). 

The frozen magnetic field is 

%rozen = 'B*{tk + At). 

This frozen field is constant in time during the whole time step. This trick was first introduced by Temam in the 
context of trilinear forms and Navier-Stokes equations for magnetic equations |26[ \T7\ , see also [TS] . 

One notices that equation b) is now a linear one, even if equation c) is still formerly non linear because it has a 
He dependence. However an explicit procedure allows to compute the solution. Indeed solutions of equation c) are 
such that 

dtm = 0. 



It means that nj and rig are frozen quantities 



ni = et 



frozen 



One has 



which yields 



dtjii u/ = J A d, d = 



frozen p> 
"'I ^frozen 

„ frozen 



nj u/ = ni ui{tk) + / V A B{s)ds A d. 



It shows that nj uj integro-differential and linear with respect to B. Plugging this form of nj u/ in equation 
b), we end up with a linear equation for B. This linear integro-differential equation is well posed under general 
assumptions. 

Once the magnetic field is computed, we can report the current J in equation c) which is now easily solved with 
the method of characteristics. It is immediate that, due to the resistive operator, the energy decreases during this 
step. Since rn is constant, the electronic energy is constant. 

The solution at the end of this stage is referred to as 

f°itk + At) and B°(ifc + At) B(tfc). 



5.3 Magnetic part, second stage 

It remains to solve 

dB 



with prescribed initial data 



8t-'^ 
dt dv 



A B) +v AB 



= 0, 



(b) 

(c), 



ntk) = f°{tk + At) and B'itk) = B°{tk + At). 



Since the magnetic field is frozen 



B = B 



frozen' 
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equation c) is greatly simplified. Once again the ionic density n/ and the electronic density rig are frozen 

ni = n/fro^^'i and = r^e^™^''''. 

A consequence is 

Cj^^frozen \ 
- ^^frozen + ^ 1 ^frozen- 

The solution of this linear equation is immediate. Therefore n/U/ is known. And finally the method of characteristics 
can be used to solve c). The energy is preserved during this second magnetic stage. 

5.4 Iterations 

The previous procedure allows us to design an approximate solution 

f{tk + At) = f'{tk + At) and B(tfe + At) = B*(tfc + At) 

one time step after the other. With this procedure the total energy decreases and we can apply our stability analysis 
for proving existence. 

A Control on moments of // 

Several type of controls on velocity moments of // are available, see [23]. Here we recall one of the most fundamental 
control in spaces based on the kinetic energy. 

Lemma 1. Let f <E .^((0,T) x n x R^) D {0,T; Ll.^{n x M.^, \v\^dxdv). Define m by ^ and n/u, by 
Then ni G {0,T] Ll^^{n)) , n/U/ e {0,T; Ll^^{n)) and we have for all t £ [0,T], 

\\niit,-)h/3<C\\f{t,;-)\\J^' ( [ [ fit,x,^)\wfdvdx)^ , (33) 



||"/(t,-)u/||5/4<C'||/(t,-,-)||oo'/'f / / f{t,X,v)\v\'d^dx) 

\Jn JR3 / 



4/5 



(34) 



Proof. We only recall the proof of the result on nj. Let R > 0. We have 

ni(t,x)^ / /(t,a;, v)dv+ / /(t,a-, v)dv 

J\v\<R "'|v|>fl 

<CR'\\f{t,;-)\\^ + ^ [ /(t,a;,v)|vpdv. 

^ "'|v|<7? 

Then by minimization over R we get 

niit,x)<C\\fit,;-)\\%'' (^|^Jit,x,v)\v\'d^r^ ' 
and after integration we obtain the claim. □ 

B Uniform lower bound on rig 

The purpose of this section is to prove several properties that we have used throughout the paper for the elliptic 
equation 

-A^Alnrie + n^, = uj, x £ Cl, 

?^ = 0, xedQ, 

ov 

with a right hand side data satisfying nj > 0, J nj = Minuiai > and nj £ L^l^ . 
We are going to prove the estimate 
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L/emma 1. Let ni G L^{0,T;L^/^ D L^ifl)) and He a strong solution to the above equation, then we have the 
two-sided control 

< A'_(||n,||5/3) <n,< A'+(ll"/||5/3), (35) 
for some continuous positive functions K±(-) with > 1 increasing, K- decreasing. 
We also recall that integration of the equation gives the electric neutrality relation 

ll"e||l = 

We can now explain why the total energy ([8]) is well defined for weak solutions and also at initial time. From 
the lower and upper bound in ([35]), we conclude that Inrie G L^(J7). Also, multiplying the equation by In rig, we 
find 

A / |VKlnne|^dx= / - rie) ln?ie < 21n/\ + j|n/||i. 

2/3 

Proof. Then, we argue in two steps. Firstly, we multiply the equation by rig and integrate by parts. The Holder 
inequality gives 

^ / ^ + / nl'^ = J njnl'^ < ||n,|U/3 ll-ell^;^, 
from which we conclude the bound 

||'^e||5/3 < \\nih/3- 

Secondly we use the elliptic regularity theory to conclude that In rig — {\nne)n G W'^'^, 1 < p < 5/3 and thus 

112 

Inrie - (lnne)j2 e L'', Vq > 1, - = -, 

q p 3 

where {(j>)n denotes the average of the L^{Q) function (j) over ft. Finally, because | > | we conclude from the 
Morrey estimates [T3] that 

II Inrie - {lnne)n\\oc < Cdlri/Hg/g). 
The result follows immediately thanks to the control of {\nnp)n through □ 

C Compactness of the magnetic field 

We have also used the Sobolev injection for Maxwell equations and we recall it in this appendix. We introduce the 
following spaces 

Hcuri(f^) = {b G L^{n)y V A b G L2(0)3}, 
Hdiv(f^) = {b G L^{^}f/ V • b G L^{nf}, 
XAr(r2) = {b G Hcuri(^^)nHdiv(f^)/ bAn = on dU}. (36) 
We recall the following result (see [9]) 

Theorem 2. Assume that the domain O is of class C^'^ . Then the space Xj^(^l) is continuously imbedded in 
Hi (17)3. 

D Kinetic averaging lemma 

We recall here one result of the theory of averaging lemmas for kinetic equations. When /(i, a;,v) is solution of 
a kinetic equation, it cannot be more regular that the initial data or the right hand-side. However, averages in 
velocity gain regularity. Recall that the macroscopic quantity is defined as 

{fm.x) = I /(t,a;,v)V'(v)dv 

where i/; is a given function in C^(R'^) (i.e. smooth with compact support), the averaging lemmas aim at proving 
compactness properties on (/'!/')■ The first version of these averaging Lemma has been established by Golse, Lions 
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Perthame and Sentis [T7] for f E ,^ ^ solution of the equation dt f + v • Vxf = S with S G ^ it was proved 

1/2 

that locally in time (fip) G iJj ^ . This version has been then be extended (by complex interpolation) for the 
framework, 1 < p < oo in |12| . Moreover, more complex versions have been proved by Di Perna, Lions and Meyer 
|12| , who treat the case where S is the fc-th derivative in velocity with a fractional derivative in x strictly less than 
one. 

We use here a version proved by Perthame and Sougadinis |24| (the equality is the exponent is due to Bouchut 
[5]) which expresses an optimal gain of regularity (a full derivative) 

Theorem 3. Let 1 < g < oo and f, g~ {gi, -.-jgd) belong to ^ .^,(R^+'^+'^) and satisfy 

§f +v.V./=|^.((I-A,.,v)^/^g). (37) 
Then {ftp) G L(^(R"'^+'^) and there exists C{p,'ip) such that 

for a positive exponent a < 5 min(i, 1 — i). 

In fact this theorem also proves strong compactness, which is the way we use it in our context. Consider a 
truncation function x — xit^^) with suppx C (0,T] x fi which is fixed in this paragraph. Define the sequence of 

functions w"^ as w'^ = yF" f^ and = Zf, with Z = ^ + v.V^cX- They satisfy in 

ot 

-(xr)+v.V,(xr) = ^.w- + z^ (38) 

We know that is bounded in H Lf ^ .^,(R^+'^+'^) for some p > 1 (here p = 30/29) and z'^ and x/^ are bounded 
inLinL,°;=,^,(Ri+3+3). 

We are going to prove the following result for an exponent q (with 1 < q < p). 

Lemma 3. Consider, after extraction, the weak limit f* in Lf ^ .^,(R^+■^+■^) of the sequence f^ weakly then 

PV- =x(r^)^p;=x(rV'), strongly in L«,^,(R1+3+3). (39) 

Proof. We first define 



= (I - A^,t,„ 



1-1/2, 



By regularizing effects the family is compact in 1^1 ^ ^{R}^^^^) : indeed g^ is given by a convolution product 
between and the fundamental solution of (I — l\x.t,v)^^^'^ ■ Moreover g^ is bounded in Lf ^ .^,(Rl+■^+■^) and then, 
by interpolation argument, it is compact in IJI^ ^(B}^^^^) for some q (with 1 < q < p). Since the family {xf^} is 
uniformly bounded well as z^, the theorem recalled above (the term z^ can be written as a v -derivative 

without loss of generality) applied to xf^ allows to claim that the following bound holds 



||,<C(g,V)C||g1l^. 
Then, there exists a subsequence f'^,w'^ and functions /*,p^,g*,w* such that 

'f^f*, w^^vf* weakly in Ll^^^{R^+^+^), 
p%^p% weakly in L«,^,(Ri+3+3), 
g' g* strongly in L^^ .j,(Ri+^+^). 
and g* = (I — Ai,_t ,u)~i/^w*. Of course we have also, passing to the weak limit, 

= ^•((i-A,,i,„)i/v) + ^r 



16 



We combine this with (|38p and we get 

+ V • v.)[xr - x/1 = I; • ((I - A..t,„)^/2(g^' - g*)) + z{r - n- 

And according to the previous theorem (the term Z{f^ — /*) can be written as a v-derivative also) once again we 
see that 

That is to say, the property holds. □ 
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